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Abstract 

We adapt our previous results for the "partition function" of the superin- 
tegrable chiral Potts model with open boundaries to obtain the corresponding 
matrix elements of e~ aH , where H is the associated hamiltonian. The spon- 
taneous magnetization Ai r can be expressed in terms of particular matrix 
elements of e~ aH S{e~^ H , where Si is a diagonal matrix. We present a con- 
jecture for these matrix elements as an m by m determinant, where m is 
proportional to the width of the lattice. The author has previously derived 
the spontaneous magnetization of the chiral Potts model by analytic means, 
but hopes that this work will facilitate a more algebraic derivation, similar to 
that of Yang for the Ising model. 
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1 Introduction 

In a previous paper [T], we obtained the partition function Zq (here referred 
to as Z p ) of the superintegrable chiral Potts model with open boundary con- 
ditions. It is a simple product of elements of two- by- two matrices, reflecting 
the fact that there is a reduced representation in which the transfer matrices 
have a direct product structure, similar to that of the Ising model. [2] 

Very recently, we have considered the problem of calculating the spon- 
taneous magnetization M. of the square lattice Ising model. [3] We used^the 
method of Yang[l] and defined M. in terms of the partition function W on 
a cylindrical lattice of L columns with fixed-spin boundary conditions on the 
upper and lower rows, with a single-spin operator Si acting on a spin located 
within the lattice. For convenience, we took the limit when the transfer matrix 
could be replaced by the exponential of an associated hamiltonian. 

The Clifford algebra technique of Kaufman [5] can still be applied to this 
system, so that W can be calculated as the square root of an L-dimensional 
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determinant. This can be further reduced to a determinant (without the 
square root) of dimension approximately L/2. 

Here we write down corresponding definitions of W for the iV-state super- 
integrable chiral Potts model, which reduces to the Ising model when N = 2. 
We conjecture in (|7.ip - (17, 3p that W is also given by a determinant of dimen- 
sion smaller than L, being a fairly immediate generalization of that for the 
Ising case. If true, this is an exact formula for finite lattices, containing three 
additional arbitrary parameters a, (3, x in addition to N, L and the labels p, q 
of the appropriate sub-spaces. It is therefore easy to test numerically, and we 
have tested it to 60 or more digits of accuracy for various small values of N, L 
(up to N + L = 10). 

If this conjecture is indeed true, then the spontaneous magnetization of 
the superintegrable chiral Potts model is given by the expression (|7.1ip below. 
This necessitates taking the limit L — > oo . As yet we have not done this, but 
we have observed numerically that (|7.1ip does indeed appear to converge to 
the known result (15. 5p . The author has previously derived (|7.1ip by analytic 
methods[6l [7] that apply in the large-lattice limit, but it would still be inter- 
esting to have an algebraic derivation that could give greater insight into the 
properties of the model on a finite lattice. 

2 Partition function 
Definition 

The model is defined on the square lattice, rotated through 45°, with M + 1 
horizontal rows, each containing L spins, as in Fig. 1. 




a a a a 



Figure 1: The square lattice C turned through 45°. 

We impose cylindrical boundary conditions, so that the last column L is 
followed by the first column 1. At each site i there is a spin <jj, taking the 
values 0, 1, . . . , N — 1. The spins in the bottom row are fixed to have value a, 
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those in the top row to have value 0. Adjacent spins a{,aj on southwest to 
northeast edges (with i below j) interact with Boltzmann weight W(<7j — aj); 
those on southeast to northwest edges with weight W(o", — 
These W, W are the Boltzmann weight functions: 

n 

W(n)= W(n + N)= p n ^[(1 - Jy)/(l - J x) , 

3=1 
n 

W(n)= W(n + N)= n~ n ]J(uj - co j x)/(l -Jy) , (2.1) 

3=1 

where oj = e 2?n / n , x and y are complex parameters, and 

= (x N - l)/(y N - 1) . (2.2) 
An important associated parameter is 

k' = ( x N -l)(y N -l)/(y N -x N ) . (2.3) 
The partition function, which depends on a, is 

z - = e n - ^) n - *■) , (2.4) 

the products being over all edges of the two types. The sum is over all values 
of all the free spins. The partition function can be written as 

Z a = ulT M u , (2.5) 

where T is the row-to-row transfer matrix, with elements 

L 

T(t,<t> = n W (^-^ + i)W(^-(rO , (2.6) 

i=i 

a being the set of all spins a%, . . . , ox in one row, and a' being the set in the 
row above. Thus T is an A^^ by A^^ matrix. The vector u a is of dimension 
N L , with entries 

(u a ) a = 1 if <7i = • • • = a L = a , 

= otherwise . (2-7) 

The sup erintegr able chiral Potts model is a special case of the more general 
solvable chiral Potts model, which satisfies the star-triangle relation. [8 J This 
ensures that two transfer matrices T, T', with different values of x, y, but the 
same value of k' , commute. 

The spin- increment matrix R 

Let R be the N L by iV^ matrix with entries 

L 

R a,a' = 116(^,^ + 1) , (2.8) 

3=1 
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where 5{a, b) = 1 if a = b (modulo N), else 5(a, b) = 0. Then pre-multiplying 
by R has the effect of increasing all spins by 1 (modulo N), hence Ru a = u a+1 
and R commutes with T: 

RT = TR . (2.9) 
For this reason it is natural to use the Fourier transform of u a : 

N-l 

Vp = N~ 1 ' 2 Y, ^ aPu a • (2-10) 

a=0 

taking p = 0, . . . , N — 1. This p replaces the Q of [SJ [I]. Then 

Rv p = u p v p . (2.11) 



If we also define 



N-l 

Z p = " pa Z a , (2.12) 

a=0 



then 



Z p = N^vtT M u 



J v — ly u p J 
N-l 

= Y 4 TM v q ■ (2-13) 

From (|2.9p , we can replace T M in the summand by R~ l T M R, and from ()2.1ip 
this is equivalent to multiplying by uj q ~ p . This in turn means the summand 
must vanish unless q = p, so 

Z p = v\T M v v . (2.14) 



The sub-space V p 

Following the observations of Albertini et aZ[lO], we showed in refs.[9j[T] that 
if one operates on v p by any product of matrices T, with different values of 
x,y but the same value of k', then all the vectors generated lie in a vector 
space V p , where p = 0, . . . , N — 1. For any vector v in V p , 

Rv = u p v . (2.15) 

We also showed that the transfer matrices satisfied a functional relation 
that determined their eigenvalues, and derived the result (12.241) for the parti- 
tion function Z p . 

If 

, \{N-l)L-p 
m = nip = integer part of 



N ' (2 ' 16) 

then there are just 2 m distinct eigenvalues. What we have not shown, but 
believe to be true, is that each such eigenvalue occurs just once, so that V p is 
of dimension 2 m . Certainly, by continuity from the case k' = 0, the largest 
eigenvalue (which is the one we most often consider) occurs just once. 
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For the case when p = and L divides by N, Au-Yang and Perk have 
recently obtained the eigenvectors explicitly. 

Two vectors v,w in different spaces V p , V p (with q ^ p) are necessarily 
orthogonal, i.e. v^-w = 0. 

Define 

N-l 

P(z N ) = z- p ^ {L+p)n {z N -l)/{z-u: n ) L . (2.17) 

n=0 

Then P(w) = P v {w) is a polynomial in w of degree m. Let its zeros be 
> • • • j w m and define 9%, . . . , 9 m 

= 0{p,l),... ,6{p,m) by 

cos% = cos[d(p,j)} = (l + ti;j)/(l -Wj) , 0<6»i<7r , (2.18) 
for j = 1, ... , m. Set 

G = (iV-l)/(^-/) , (2.19) 

5 = JVCl-x-^/a-x-") , (2.20) 
define the two-by-two matrices 

1 - x~ N 

F(x,y,9) = — - [GI 2 + (1 -k' cos 9)S -k' sin0 C] , (2.22) 

I2 being the identity matrix, and set 

D(cos9) = ( 1 )-F(x,y,6) M - ^ (2.23) 

then in [T] we find that 

Z p = g LM x~ Mp D{cos9i)D{cos9 2 )---D(cos9 m ) . (2.24) 

3 The hamiltonian limit 

Take 

H = e" 2e (3.1) 
and consider the limit when e — > 0. Then to first order in e 

x = y = 1 + 2k' e , (3.2) 

W(n) = l-2ne , W(n) = 2fe'e/(l - w~ n ) (3.3) 
for < n < iV, while W(0) = W(0) = 1. Noting that 

7V ~ 1 

JV - 1 - 2j = 2 V (3.4) 

71=1 
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for < j < N, it follows that 

T = [1- (JV-l)Le] I-eU , (3.5) 
where I is the identity matrix and 

L N-l 

H = -2Y / Y, {Z?Zj?i + k'X])/(l - . (3.6) 

j = l 71=1 

This is the hamiltonian associated with the transfer matrix T. Since all trans- 
fer matrices with the same value of k' commute, they also commute with H. 



Here Zj , Xj are the N L by N L matrices of |10j , with elements 



{Z 3 )^,=^ Y[5(a m ,a' m ) , (3.7) 

m=l 

L * 

(Xi)a*> = °j + 1) II S ( a n, <) , (3.8) 

n=l 

the * on the last product indicating that that it excludes the case n = j. 

The hamiltonian 7i is known to have very special properties. In particular 
Au-Yang and Perk showed that it satisfies the "Onsager algebra" . |12j 

Still working to first order in e, we obtain 

g = l + (N-l)k'e , 

{1 ~f k , N)G = l-N(l + k>)e , 
F(x, y, 0) = [l- N(l + k')e]I 2 + Ne [(1 - k' cos 9)S - k' sin 9 C] . 
Now we take 

e = a/M , (3.9) 
and let M — > oo, keeping a fixed. Then 

F(x, y, 6) M -» exp{iVa[-(l + k')I 2 + (1 - k' cos 9)S - k! sin 9 C}} (3.10) 



and from (|3.5p . 

T M _^ e _(jv-i)ia exp (_ aW ) . (3.H) 
From (pl23j) and flQlD , it follows that 

v\ exp(-aH) v p = e^ a S(cos6»i) ■■■D(cos0 m ) , (3.12) 

where 

/i = /ip = 2/c / p + (l + A; / )(^A r -^ + ^) , (3-13) 
D(cos9) = ( 1 ) • exp[-aF(0)]- ( J^J (3.14) 

and F{9) is the two-by-two matrix 

F(9) = -N(l - k' cos 9)S + Nk' sin 6* C . (3.15) 
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The two-by-two exponential 



We can calculate the exponential in (|3.14p of the two-by-two matrix —aF(9) 
in the obvious way, by diagonalizing it, exponentiating, and then returning to 
the original basis. If we define 



A = A(0) = (1 - 2k! cos 6 + k' 2 ) 1 / 2 , 
u p (a,6) = cosh(NaX) + sinh(iVaA) 



v p (a,9) 



k' sin 9 
A 



A 



sinh(iVaA) 



(3.16) 



(3.17) 



w p (a,9) = cosh(NaX) — sinh(iVaA) 

A 



then 



Hence 



u p (a,9) v p (a,9) 
v p (a,9) w p (a,9) 



exp[-aF(9)] = 

D(cos9) = u p (a,9) 



and (|3.12p becomes 

u^exp(— aH)v p = e~ flpa u p (a,9i) ■ ■ ■ u p (a,9 m ) 

4 Reduced representation of 7i 



(3.18) 
(3.19) 
(3.20) 



We consider some basis of the 2 m -dimensional vector space V p and label the 
vectors by s = {s%, . . . , s m }, where each Si takes the values 1 or —1. We 
can think of the Si as "Ising spins". Thus there are 2 m vectors v s = v% = 
v(si, . . . , v m ), each of dimension N L . 

In pQ we showed that we can choose the vectors v s so that v p above is 



and 



Hv(si,. . . ,s m; 



v p = u(l, !,...,!) , 



fi- N^2(l- k' cos 6 >j)sj 
i=i 



(4.1) 



v(su ...,s m ) + 



l , . . . , Sj , . . . , s r 

3=1 

Defining 2 m by 2 m matrices Sj , Cj by 

m 



(4.2) 



n=l 



(Cj) 



6 ( S ji- S 'j)n* 6 ( S n,s' n ) 



(4.3) 
(4.4) 



n=l 
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where again the * means that the term n = j is excluded from the product, 
we see that with respect to this basis the hamiltonian 7i is now 



H = v P - Nj2[(l-k' cos 9 j)Sj -k'smOjCj] , (4.5) 
3=1 

which is equation (2.20) of pQ. This is consistent with our result fj3. 121) above. 
From (pn3|) . (|4"3|) . H is linear in k! . Set 

H = H + k'H 1 , (4.6) 

Hq,Hi being independent of k' , and define 

in 

k(s) = £(l-s;)/2 ■ ( 4 - 7 ) 

3=1 

then k(s) takes the integer values 0,1, ... ,m. If we order the rows and columns 
of H with increasing values of n(s), then Hq is diagonal and Hi is block tri- 
diagonal, with non-zero entries only when \n(s) — k(s')\ < 1. 

From (14. 5p . H is a direct sum of m two-by-two matrices. Similarly, if we 
define the two-by-two matrix 

u = ( u p (a,ej) v p (a,0j) \ 

then exp(— aH) is the direct product 

exp(-aif) = e~ ra U x ® U 2 <8 • • • 8> i/ m . (4.9) 

Let |0) be the 2 m -dimensional vector whose elements s are zero except for 
the element si = 1, s 2 = 1, ■ ■ ■ , s m = 1, which is unity, i.e. 

|o>-(jVj)« ■■■•(;) ■ W.10) 



«}e- aW «„ = (0|e- Qi/ |0) (4.11) 



This is the representative of the A^ L -dimensional vector v p . If (0| is the trans- 
pose of |0), then 

and equation (|3.20p follows immediately. 

The derivation of [9j [1] does not exclude the possibility that the basis 
vectors v s depend on the parameter k' . However, all studies for small N,L 
agree with the hypothesis that they are (or at least can be chosen to be) 
independent of k'. This is consistent with the fact that both 7i and H are 
linear in k' . 
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5 The spontaneous magnetization. 



Consider the lattice of Figure Q] and take a = 0, so all upper and lower 
boundary spins are fixed to be zero. Let £ be the spin on a site deep inside 
the lattice. Then in the usual way we can define the order parameters of the 
chiral Potts model as 

Mr = (u/" C ) (5.1) 

for r = 1, . . . , N — 1. Here (/(C)) denotes the usual statistical mechanical 
average 

(/(C)) = Z, 1 ]T /(C) II W{ai - aj) J] W(a t - *j) (5.2) 

for any function /. We take the limit when the lattice is infinitely large, so 
L,M — > oo, and C is infinitely far from the boundaries. 

The W, W products are unchanged by incrementing all spins by one, so if 
we imposed toroidal boundary conditions, then it would be true that 

(/(C + 1)) = (/(C)> (5.3) 

and this would imply that M r = uj t M r - Hence for (mod N) we would 

necessarily have M r = 0. 

At high temperatures {k' > 1), this is true also for our fixed-spin boundary 
conditions when we take the large-lattice limit. However, at lower temper- 
aturers (0 < k' < 1) the system has ferromagnetic long-range order and 
"remembers" the boundary conditions even in the limit of C deep inside a 
large lattice, and 

< Mr < 1 • (5.4) 
In fact we know M r - In 1989 Albertini et al |10j conjectured that 

M r = (1 - k >Y N - r)l2N2 (5.5) 

and the author was able to derive this formula in 2005 [HI \7\. The method used 
was analytic, depending on the star-triangle relation, functional relations and 
analyticity properties. 

When N = 2 the chiral Potts model (both superintegrable and general) 
reduces to the Ising model, whose partition function was obtained by Onsager 
in 1944. [13j Onsager announced at a conference in Florence in 1949 that he 
and Kaufman had solved the spontaneous magnetization and obtained M\ = 
(1 — A; /2 ) 1 ^ 8 ,[14J but the first published derivation of that result was given by 
Yang in 1952.(4] 

Onsager and Yang's methods were much more algebraic, determining the 
eigenvalues of the transfer matrix T, and certain elements of the eigenvectors. 
It would be interesting to obtain a derivation of M r that parallels Yang's. 
The object of this paper is to suggest how one may make progress in that 
direction. 

We introduce the N L by iV^ diagonal matrix S r with elements 

L 
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Note that, for all integers p and r, 

SrVp+r = V p , V p ^S r = V^ +r . (5-7) 

Because of the cylindrical boundary conditions, we can take the spin £ 
to be in any column, so we choose it to be in column 1. Then (|5.ip can be 
written 

M r = W/Z , (5.8) 

where 

W = ulT j S r T M - j u , Z = uIt m u , (5.9) 

j being the number of rows below £. 
Prom OTTO]) and ([27L2D . 

N-l N-l 
W = N- 1 Y J vlTiS r T M -i Vq , Z = N- 1 J2 4T M v p . (5.10) 

p,q=0 p=0 

Since R commutes with T and 

RS r = oj"'S r R , (5.11) 

it follows from (|2.15p that the first summand in (|5.10p vanishes unless q = p+r, 
so 

Af-l 

W = N~ 1 J2 vlTiS T T M -j V p +r , (5.12) 

p=0 

interpreting p + r asp + r to modulo N. 

For < k' < 1 and L is large, the N largest eigenvalues of T are asymp- 
totically degenerate, their ratios being of the form 1 + 0(e~ Lu ), v being a 
measure of the interfacial tension. However, there is one and only one of 
these eigenvalues in each of the vector spaces V p , for p = 0, . . . , N — 1. 

Since T and TL commute and TC is hermitian, the eigenvectors ip p corre- 
sponding to these eigenvalues are unitary, so 

V^q = 5p,q ■ (5-13) 



Asymptotic degeneracy 

In each sub-space V p there is single largest eigenvalue A p of the transfer matrix 
T and these eigenvalues are asymptotically degenerate, in the sense that for 
large L there is a common value A such that 

A- X A P = l + 0(e~ Ls r) , (5.14) 

i.e. the ratios of the A p approach unity exponentially rapidly. 

This can be seen by considering the series expansion of the eigenvector ip q 
in powers of k' . Since T, TL commute, we can look at the eigenvectors of TL, 
corresponding to the most negative (ground state) eigenvalue. 

When k' = 0, TL = TL , where 

L N-l 
j=l n =l 
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This is diagonal, with minimum eigenvalue — 2L, when all the L spins are 
equal. Thus from (|2.7p . no, . . . , Ujv-i are ground state eigenvectors. 

We can start from one of these eigenvectors and use standard linear pertur- 
bation theory to develop a series expansion for the eigenvector of TC, starting 
from the initial eigenvector u a . This entails changing successively more of the 
spins from value a to some other value. It will work until all of the spins 
are changed, when for the first time we come to another of the eigenvectors 
of TLq. At that stage, and only at that stage, one would have to resolve the 
degeneracy of the initial eigenvalues. This means that naive perturbation 
theory works to order k' L . The calculation only depends on a in so far as it 
involves the differences (mod N) of the L spins from a. Thus to this order 
the eigenvalue is independent of the initial choice of a. This is true also of the 
eigenvalues of T, so A p = A, A being the common eigenvalue, in agreement 
with (157HD . 

Also, if ip' a is this near-eigenvector, then 

tpa ^b = &a,b , (5.16) 

where £ is independent of a and b, and to this order the actual eigenvectors 
are 

N-l 

VV = N- 1 ' 2 J2 "~ aP < • (5-17) 

a=0 

it follows that, for all p, 

4 v p = t • ( 5 - 18 ) 

In the limit of j,M — j,L large we can replace T J in (|5.10p . (|5.12p by 
ippA^tp^ (with the appropriate value of p), and T M_J by ^pA^ - - 7 '^, giving 

vlTiS r T M -i Vp+r = e MA eZ4Srip P+ r , vlT M v p = e MA Ct (5.19) 
£* being the complex conjugate of £ and 

ip^S r ip p+r = independent of p . (5.20) 
Thus W, Zq are the two expressions in (|5. 19|) , respectively, and 

M r = ^Srlpp+r . (5.21) 



Expressions in terms of 7i 

Rather than continue to work with the transfer matrix T, we find it convenient 

to instead use the negative exponential of the hamiltonian and to replace 
T j^ T M-j in (py^j) by e -«W e -/JW and T M in dgjQj) by e -aH ( with a different 

a), making them 

N-l N-l 

W = N- 1 ^ W m , Z = N- 1 Z g , (5.22) 

where now, setting q = p + r, 

W p ,g = W p , g {a,f3,x) = vle- aH e- pJ S r e-P n v q (5.23) 
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Z p = Z p (a) = V y aH v p (5.24) 

and 

x = e~ 2N P . (5.25) 

We have introduced the matrix factor e~ p ^ immediately pre-multiplying 
S r in (15T23J) . Here 

J = H + L{N-1)I (5.26) 

is a diagonal matrix whose entries are 0, 2N, 4N, . . . , 2N[(N — 1)L/N\. Hence 
W Pj g(a,P, x) is a polynomial in x of degree [(N — 1)L/N]. This naturally 
manifests itself in the following working and provides a useful check against 
errors. 

We can think of these Z p , W Ptq as hamiltonian partition functions. They 
are rather simpler than the original partition functions to work with. 
When p —* + oo, then x — > and e~ p ^ — > v p v p \ so, using (|5.7p . 

= Z p (a)Z q (0) , (5.27) 

Wj,,,(a,0,!r) = Z p (a) , W p>q (0,p,x) = Z q {f3) . (5.28) 

These relations also provide useful checks on our subsequent calculations. 

Because Tt, T commute, they have the same ground-state eigenvectors ip p . 
In the limit when p = 0, and a, (3, L — > oo, we obtain 

W m (a,(3,l) =e-^ A t^4Sr^ q , 

Z p (a) = e" aA re • (5.29) 

So from H52QJ), (I53TD . 

^ r = lim (5.30) 

a,/3,L-.oo (Z p (2a)Z,(2/3))V2 

for any p, q such that < p, q < N and q = p + r, mod iV. 
From (I330|) and (I5T24T) . 

Z p (a) = e-^ a u p (a,0 1 )---u ¥ (a,0 m ) . (5.31) 

It remains to calculate W p q (a,0,x). We have not done this, but the rest 
of this paper is concerned with presenting a conjecture for it as a determinant 
of dimension not greater than (iV — 1)L/N. This expression agrees with the 
known N = 2 result for the Ising model, and indeed is a fairly immediate 
generalization of that result. It has the properties (I5.27|) . (15.281) . and has 
been extensively tested numerically for small values of N,L. 
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Expressions in terms of H 

First we remark that if v E V p + r and v' = S r v, then from (|2.15p . Rv' = uj p v', 
so v' is a candidate for the sub-space V p . However, in general it does not lie 
within this sub-space. Even so, we can define a matrix S^ ed of dimension m p 
by m p+r by 

(S^y = (v^SrV^ . (5.32) 

These elements depend on N,L,p,r. They are of course independent of a 
and (3. From our remarks at the end of section 4 that we expect the v s to be 
independent of k', the same must be true of the elements of S^ cd . 
We can then write (|5.23p as 

W M = (0\e- aH e-» J S: ed e-P H '\0) , (5.33) 
Z p = (0\e- aH \0) , 

where H' is the H of (|4.5p . (|4.6p but with p replaced by p + r, and 

m 

J = H + L(N - 1)1 = Nj^i 1 ~ s j) ( 5 - 34 ) 

3=1 

is the diagonal matrix with elements 2Nk(s) in position (s,s). 
Let 

u p ( 1, a, 9) = u p (a, 0) , u p (—l,a,9) = v p (a,0) . (5.35) 

and set 

q=p + r , m' = m q , // = fj, q , 9' j =0(q,j) . (5.36) 
Then we can write these equations more explicitly as 

Wp, q (a,f3,x) = e' 011 *- 13 * 1 ' ^Up(s 1 ,a,9i)---u p (s m ,a,9 m ) X 

s,s' 

x^Hs: cd ) s ^u q (s[,(3,0 > 1 )---u q (s' m ,P,9' m/ ) , (5.37) 

and 

Z p (a)=e- a ^u p (a,9 1 )---u p (a,9 m ) . (5.38) 

The non-zero elements (s, s') of the 2 m by 2 m ' matrix S^ ed satisfy k(s) = k(s'). 
If we also order the rows and columns of S^ ed in increasing value of k(s), then 
this matrix is block-diagonal. 

We do not have a direct derivation of S^ ed , though of course it can be 
calculated numerically for small values of N, L from (|5.32p . In principle it can 
be calculated from our conjecture (|7,2p below. If s is the m by m' diagonal 
blocks of S'£ ed in the block k(s) = n(s') = 1, h is the corresponding m by m 
block of Hi, and h' the m' by m' block of H[, then this conjecture implies 
that the double commutator h-h-s — 2 h-s-h' + s-h'-h' is of rank one. This was 
a key initial encouraging observation in our search for the expression (|7,2p . 
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6 The orthogonal matrix B 

Before stating our conjecture, we define an m by ml real orthogonal matrix 
B = B pq whose elements involve the 9\ , . . . , 9 m defined by fl2. lTf) , (|2.18|) , as 
well as the 9[, . . . , 9' m , defined similarly, but with p replaced by q and m by 
ml . We must have p ^ q. 

We define B = B Ptq to be the matrix with elements 

Bi,j = f {p,q,i)f{q,p, j)/ (cos 9i - cos 6'j) , (6.1) 

where we choose the functions f(p,q,i),f(q,p,j) to ensure that 

B T B = I if m > ml , BB T = I if m < rri , (6.2) 

/ again being the identity matrix, of dimension min(m, ml). 

The case p < q 

From (|2.16p , if p < q, then m > ml and we want B T B = I. From (|fTLj) . 

_ f{q,p,i)f{q,pJ) y> f f(p,q,n) 2 f(p,q,n) 2 \ , g3 , 

cos 9'a — cos 9'i , I cos 9' 4 — cos 6 n cos 9'- — cos 9 n \ 

for i ^ j. 

We want the RHS of (16. 3p to vanish for i 7^ j. Consider the functions 



P P (c) = U(c- cos 9 t ) = N- L (c + l) m p(£ [ ) , (6.4) 



\c+l J 
i=l 



' ■ /" — nr\c H 



C — COS n 



n=l 

The first is a known function, given by (|2.17p and f|2. 18|) . the second is of the 
form IZp(c) / P p (c) , 1Z p (c) being a polynomial of degree m — 1. We want there 
to exist constants 7, 7' (dependent on p, q) such that 

T(c) = j' + 7 P q (c)/P p (c) , (6.6) 

since then .F(cos^) = IF(cos9j) = 7' and the RHS of (|6.3p vanishes. This 
implies that 

K p {c)=iP p (c) + 1 P q {c) . (6.7) 

From (|2.16p . m and ml = m q differ by at most one, so ml + 1 > m > ml . 
Whether m = ml or m = ml + 1, we can always choose 7' to ensure that the 
RHS of (|6.7p is a polynomial of degree m — 1. Then the equation defines 1Z p (c) 
(to within the factor 7) and the parameters f(p,q,n). 

From (I6.5p . /(p, g, n) 2 is the residue of ^"(c) at the pole c = cos 9 n , so from 

f(p,q,n) 2 = 7 Pg(cos6> n )/A p (cos6> n ) , (6.8) 



14 



where 

A P (c) = ^P p (c) . (6.9) 

For given p, q, this determines f(p, q, i) to within a factor independent of 
i (but possibly dependent on p and q). To determine this factor we need to 
consider the case when i = j in the first of the equations (|6.3p . which gives 



f( P ,q,i) 2 G(cos6i) = l , (6.10) 

where 

m 

G(c) = £ /(p, ? , n) 2 /(c - cos 9 n f . (6.11) 



n=l 

From the equations above, 



G(c) = ~F(c) = -7^- ^4 • (6.12) 



Since P q (cos9j) = 0, this gives 

G(cos^) 7A g (cos^) 

The parameter 7 is at our disposal. We observe numerically that for small 
values of n and L we can ensure that f(p, q, n) 2 , f(q, p, i) 2 are real and positive 
by choosing 

7 = 1. (6.14) 

We can then take f(p, q, n), f(q,p, i) to be positive, for all n, i. The matrix B is 
then defined by (|6.ip , (|6.8|) , (|6.13|) , (|6.14|) . It is real and has the orthogonality 
property B T B = I. If m = m! this implies BB T = I. 

The case p > q 

We can combine (|6.8p . (|6.13p into a single formula by defining 

e(p,g) = l ifp<9 , e(p,?)= -1 ifp>g • (6.15) 

Then both equations are contained in 

r -n 1/2 

f(p,Q,i) = [e(p,q) P q (cos9i)/A p (cosdi)\ , (6.16) 

for p ^ q. 

We can now extend the formula (16. ID to all p 7^ g. It is readily observed 
that 

S,, p = -P^ g . (6.17) 

We have just established that B T B = I if p < q. It follows that BB T = I if 
p > q (which implies m < m'). This is the desired orthogonality property. 

We remark that we have only conjectured (based on numerical calcula- 
tions) that the RHS of (|6.16p is real and can be chosen positive. If this were 
to fail the above formulae would still apply, but B pq would be a complex 
orthogonal matrix. 
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The matrix E 

We shall also need the mbym diagonal matrix E pq , with entries 

[ E P,q]i,j = e 0>9^) s i,j > ( 6 - 18 ) 
where the function e(p, q, i) is defined as follows, for < p,q < N: 

e {Pi Qi i) = sin#j if p < q and m > m 

= tan(0j/2) if p < q and m = m! 

= l/sin#j if p> q and m < m! (6.19) 

= cot(0j/2) if p > q and m = m' . 

Since m — 1 < m! < m if p < q, and m + 1 > m' > m if p > q, these equations 
cover all cases; 9i = 6(p, i) is again as defined in (12, 18ft , The function e(q,p, i) 
is defined similarly, but with p, q interchanged and 6i replaced by 9^ = 9(q, i). 



7 The conjecture for W 

We return to considering the W Ptq of equations f)5.23j) , f)5.33j) and (|5.37p . Based 
on the calculation for the Ising model, [3} eq.7.9] we conjecture that 

W p , q (a,f3,x) = Z p (a)Z q (J3)V M (a,P) , (7.1) 

where T> p ^ q (a,(3) is the mbym determinant 

T> P;q (a,f3) = det{I m -xX p {a)E P:q B Piq X q {(3)E g , p B g:P \ (7.2) 

or equivalently the m' by m' determinant 

V P)q {a,l3) = det[I m i - xX q (P)E q>p B qiP X p (a)E P)q B Pjq ] . (7.3) 

Again I m is the identity matrix, of dimension m and X p {a) is the diagonal m 
by m matrix whose entry in position is 

[X P («)kj = V ^jz¥)^ (7-4) 
u p {a, Vj) 

Note from (|6.ip that each function f(p,q,i),f(q,p,j) occurs twice (i.e. as 
its square) in (17.21) and (|7.3p , so the choice of the square roots in (|6.16p is in 
fact irrelevant. 

From (|3T7ft and (|fT23|h v p (a, 6) = and Z p {a) = 1 if a = 0, so flH]) does 
indeed have the properties (|5.27p . (|5.28p . It is a fairly immediate generaliza- 
tion of eqn. (7.7) of [3] and has been tested to high numerical accuracy (60 
digits or more) for arbitrary k' , a, (3 and all N, L,p, q such that 2 < N, 3 < L, 
N + L < 10. We conjecture that it is true for all N, L,p, q, x, a, f3. 
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Consequences 

Define 

a pJ = {1 - k'e 1 ^} 1 ' 2 , b Ptj = {l-k'e- ld i} 1/2 , (7.5) 

where 9%, . . . , 9 m are given by (|2.18p . They depend on p. Again the function 
m p is defined by (|2.16p for < p < N, and m = m p , m! = m q . 
Then from (I3.16|) . 

\ j = A(%) = (1 - 2k' cos Oj + k' 2 ) l/2 = a p>j b Pij , 

so from (|OID and pTTTD . 

lim 4^ = TT ^±M! . (7 . 6) 
Z p (2a) ^ 4a Pij 6 PJ 

Also define quantities x p j, not to be confused with the x of (|5.25|) . by 

*Vi = lim ^44 = x ~ k ' S ^ 3 a ■ (7-7) 



a— >oo 



Then 



and 



b 



*pj = > j^—r 1 (7-8) 



i=i 



Let 



A p = lim A p (a) , (7.10) 



a— >oo 



so it is the diagonal matrix with diagonal elements x p j. Taking the limits 
a, (3 — > +oo and setting x = 1, it follows from (|5.30p . (|7.ip that if g = p + r 
to modulo A, then 

det(/ m XpEp q Bp qXqEq p B„ p) 

M r = lim 9 7.11) 

i-oo {det(/ m + X 2 ) det^+A 2 )} 1 /? 

where < p, q < N and < r < N. 

We have not been able to evaluate the RHS of (|7.11|) analytically. Even 
for the A = 2 Ising case discussed in [3] , we do not directly evaluate (|7.1ip , 
but rather the expression in terms of square roots of L by L determinants 
that leads in that case to (|7.11|) 3 

We have conducted numerical experiments for various values of N,p, q and 
k', and observed that as L — > oo the expression on the RHS of (|7.1ip does 
indeed approach the known result (|5.5p . the error for finite L being of the 
order of k' L or smaller. 



1 We do this by writing as the determinant of a Toeplitz matrix and using Szego's theorem. 
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8 Summary 



We have defined the hamiltonian partition functions W p ? (a, j3, x), Z p (a) by 
(|5.23j) . (|5.24p and shown that the spontaneous magnetization A4 r of the su- 
perintegrable chiral Potts model is given by (|5.30p . For the general solvable 
chiral Potts model, M. r is independent of the rapidities[7j p. 7]. The superinte- 
grable model is obtained from the general by a special choice of the rapidities 
{k! being the same), [U p. 5] so M. r is the same for bothH 

By taking the hamiltonian limit of the results of [I], we show that Z p {a) 
is given by (I5.38p . We then conjecture that W p<q {a, f3,x) is given in terms an 
m by m determinant by (|7.ip . This is a natural generalization of the known 
result for the special case N = 2, i.e. the Ising model. [3] 

If this is true (and all the numerical evidence suggests that it is) this is 
a huge simplification, reducing the problem from exponential complexity to 
comparitively small polynomial complexity. Even so, we have not been able 
to make the final step and to obtain M. r from (IT. 1 1H . We already knowjBJ [7] 
that M. T is given by (|5.5p . but it would be interesting to obtain it by this 
more algebraic route. The matrices B pq and (for m > m') T> Piq (a, P)B Ptq are 
Pick matrices. |15j 

So there remain two things to do: to prove the conjecture (|7.ip and to 
evaluate the limit (|7.1ip . The first is an algebraic problem, the second an 
analytic one. The fact that (|7.ip contains the additional parameters a, (3, x 
should be helpful in establishing it. 
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